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Model study of the sign problem in the mean-field approximation 
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We argue the sign problem of the fermion determinant at finite density. It is unavoidable not only 
in Monte-Carlo simulations on the lattice but in the mean-field approximation as well. A simple 
model deriving from Quantum Chromodynamics (QCD) in the double limit of large quark mass 
and large quark chemical potential exemplifies how the sign problem arises in the Polyakov loop 
dynamics at finite temperature and density. In the color SU(2) case our mean-field estimate is in 
excellent agreement with the lattice simulation. We combine the mean-field approximation with 
a simple phase reweighting technique to circumvent the complex action encountered in the color 
SU(3) case. We also investigate the mean-field free energy, from the saddle-point of which we can 
estimate the expectation value of the Polyakov loop. 

PACS numbers: ll.10.Wx, 11.15.Ha, 11.15.Tk, 12.38.Aw 



I. INTRODUCTION 

Quantum Chromodynamics (QCD) exhibits various 
states of matter depending on the environment [l[: 
A strongly coupled quark-gluon plasma (sQGP) has 
been discovered at the Relativistic Heavy Icon Collider 
(RHIC) and various color superconductors presum- 
ably exist in the cores of the compact stellar object 
These are widely-known extreme states of equilibrated 
QCD matter at high temperature and at high (baryon or 
quark) density, respectively. The study on sQGP prop- 
erties at high temperature has been led by the Monte- 
Carlo simulation of QCD on the lattice [4(. The lattice 
QCD results provide us with fundamental information 
such as the phase transition temperature T c Q , equation 
of state [6[, susceptibility j7J, mesonic correlation above 
T c @, etc. 

In contrast, at finite density, the lattice technique is 
not quite successful so far; it is hindered by the notori- 
ous problem that the fermion determinant in the presence 
of nonzero quark chemical potential fi q is not positive 
semidefinite (i.e. not nonnegative) . This problem is com- 
monly refereed to as the fermion sign problem [9|. The 
Monte-Carlo simulation based on importance sampling 
requires positive semidefinite probability for each gauge 
configuration. In the presence of /i q ^ 0, however, the 
probability is no longer a well-defined quantity due to a 
negative fermion determinant arising. There have been 
several techniques proposed to handle the sign problem, 
e.g. the reweighting method [lfj, the Taylor expansion 
method [ll[, the analytical continuation from an imag- 
inary chemical potential [Hj]. Although these proposals 
have achieved partial successes when fi q /T is small, no 
prescription applicable at high density has been estab- 
lished yet. 

This work is aimed to point out that the sign prob- 
lem is relevant even in the mean-field treatment of gauge 
fields. One can intuitively understand it in the following 
way; in the mean-field approximation the partition func- 
tion is estimated from the most dominant contribution 
of particular configurations that have the largest proba- 



bility or the smallest free energy. When the probability 
for each gauge configuration is ill-defined due to the sign 
problem, hence, the mean-field free energy may well be 
problematic. 

In the present work we shall focus on a specific man- 
ifestation of the sign problem appearing in the simplest 
model at finite temperature and density. Because QCD is 
a highly nontrivial theory, applicability of the mean-field 
approximation is quite limited. We choose a finite tem- 
perature model because the mean-field description pre- 
sumably works well to sketch the hot QCD phase transi- 
tion. The Polyakov loop plays an essential role as an or- 
der parameter there [l3| . The dynamics of the Polyakov 
loop was closely examined some years ago both in the 
lattice simulation [l4j and in the mean-field approxima- 
tion [IE) 13- Recently, in addition, the Polyakov loop dy- 
namics near T c is specifically paid attention [13] • There 
is also an interesting observation that the entanglement 
between the chiral and Polyakov loop dynamics turns out 
to be indispensable to understand the nature of the QCD 
phase transitions [H, QjJ H(| • 

We have already known several indications from the 
mean-field studies about the Polyakov loop behavior at 
[i q 7^ 0. In the best of our knowledge, the effective poten- 
tial of the Polyakov loop at finite temperature and den- 
sity was first derived in Ref. [2l[ in perturbative QCD 
in the one-loop order. (See also Ref. [22j for the po- 
tential with an imaginary chemical potential.) It is ob- 
vious from Eq. (5) in Ref. [2lJ that the Polyakov loop 
variable is augmented to a complex valued variable when 
fj,q =/= 0, and thus the effective potential turns complex. It 
is quite nontrivial how to derive meaningful information 
from such a complex effective potential; the expectation 
value of the Polyakov loop cannot be fixed directly by a 
minimum of the complex effective potential. The same 
problem occurs also in the chiral effective models with 
the Polyakov loop coupling which was first formulated 
by one of the present authors in Ref. [2(| and has been 
investigated extensively [H, |24| . though the sign prob- 
lem has been almost overlooked. Here, we would refer 
to a closely related work, Ref. (2f| , in which the authors 
mentioned on the sign problem in a matrix model of the 
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Polyakov loop. It should be noted that the mean-field 



is not complex unlike 
so that one can deter- 



free energy as seen in Refs. [24, 
the effective potential in Ref. 
mine the expectation value of the Polyakov loop without 
ambiguity. We will go into details of this issue later. Our 
findings are all consistent with what has been foreseen in 
Ref. Ig. 

The concept of this work lies not in solving the sign 
problem but in observing the sign problem in a manage- 
able way, so to speak, in a clean environment. We also 
demonstrate that a technique which we call the phase 
reweighting method works fine as a practical prescription 
in the similar but not the same spirit as in the finite 
density lattice simulation. We will employ the dense- 
heavy model obtained from QCD in the double limit of 
large quark mass and large quar k chemical potential [26| 
(see also discussions in Ref. [27|). The reasons we adopt 
the dense-heavy model are as follows: First, the fermion 
determinant is exactly calculable as a function of the 
gauge field. Second, neither the chiral condensate nor 
the diquark condensate is involved in the dynamics ow- 
ing to lar ge q uark mass. Third, lattice data is available 
from Ref. [26| . From these three reasons the dense-heavy 
model is considered to be an appropriate implement for 
our purpose to scrutinize the sign problem within the 
mean-field approximation. 

This paper is organized as follows. In Sec. |TT]we make 
a brief overview on the sign problem. We formulate the 
model in Sec. IIII Al and the mean-field approximation in 
Sec. IIIIBI Then, in Sec. IIII CI we examine the color 
SU(2) case first which is free from the sign problem in 
order to make sure if the mean-field approximation is 
reasonable. We next proceed to the SU(3) calculation. 
In Sec. IIII Dl we explain the phase reweighting method to 
circumvent the SU(3) complex fermion determinant and 
then discuss the validity of our method by viewing the 
mean- filed effective potential in Sec. IIII El Section IIVI is 
devoted to the summary. 



II. SIGN PROBLEM 

Here are general discussions on the sign problem of 
the fermion determinant at finite density. Readers who 
are already familiar with it can skip to the model study 
starting from Sec. IIIII We will later demonstrate how 
the dense-heavy model concretely embodies the general 
features mentioned in this section. 

The fermion determinant in Euclidean space-time with 
a quark mass m q and a quark chemical potential fi q takes 
the form of 



det M(ji q ) = det^-D** + 74/Zg + m q ] , 



(1) 



where = — igA^ is the covariant deriva- 

tive. It is a well-known argument that detyV((/i g ) = 
det 75^(^)75 = {det M(— (i q )}* and thus the fermion 
determinant is real in the zero density (fi q = 0) case. 
Alternatively, we can explicitly look into the eigenvalue 



spectrum of the Dirac operator to confirm that the de- 
terminant is positive semidcfinite as follows: For fi q = 0, 
when ip n is an eigenstate of 7 M -D M , the eigenvalue A n is 
pure imaginary because 7^-0^ is anti-Hermitian where 
7 M 's are Hermitian in our convention. Since the mass 
term is simply proportional to unity, ip n is also an eigen- 
state of 7 M -D M + m q with the eigenvalue A„ + m q . One 
can show that 75?/>„ is an eigenstate of j^D^ + m q hav- 
ing the eigenvalue — A„ + m q — (X n + m q )* from the 
property 757 M -D M 7s = —j^D^. Therefore, the determi- 
nant consists of a pair of A ra + m q and (A„ + m q )* which 
makes the whole determinant real and nonnegative as 
|A„+m 9 | 2 >0. 

In the presence of the chemical potential term the 
fermion determinant is not necessarily positive semidefi- 
nite. When ip n is an eigenstate of 7^-0^ + 74M<? with the 
eigenvalue A„, in the same was as above, one can show 
that 75^ is also an eigenstate of j^D^ + 74/Xq with the 
eigenvalue — A„ which is not A* because 7 M -D M + 74^ 9 is 
no longer anti-Hermitian. It is apparent from this expla- 
nation that, if fi q is pure imaginary [l2j], 7 M £> M + 74/-% 
is anti-Hermitian so that the fermion determinant turns 
nonnegative. 

If a considered theory has degenerate fermions associ- 
ated with an internal symmetry under the transformation 
T and the chemical potential \x q is replaced by a matrix 
(j, in internal space transforming like 



T- 1 nT 



(2) 



then detM(fi) = {det A4(/j,)}*, so that the determinant 
becomes real. 

A well-known realization of this comes from the isospin 
degrees of freedom in flavor (u,d) space. The isospin 
chemical potential fi a 73 and T — iti (or T = n) 
certainly satisfies Eq. |(5J) where r's are the Pauli matrices 
in (u,d) space |28j . 

Another realization is the color SU(2) case [2!| that 
is relevant to our model study as we will argue later. 
The C-transformation changes the quark chemical po- 
tential fig as in Eq. @. In general cases, however, the 
C-transformation does not give an escape from the sign 
problem because it also swaps color and anticolor. That 
is, 



C- 1 75 7 / ^7 5 C= b^-igW)]' 



(3) 



with (A^) c = -(A^)*. Special for N c = 2 is that anti- 
color is not distinguishable from color since two doublets 
can make a singlet. Actually <72(A fl ) c 02 = A* 1 where 
it's are the Pauli matrices in color space. Therefore the 
SU(2) case is free from the sign problem. 

We should note that the fermion determinant being 
complex is not necessarily harmful on its own. Rather, 
it is important whether the real or imaginary part of 
the fermion determinant is positive semidcfinite or not. 
Even though the determinant evaluated for a certain A M 
is a complex number, the functional integral over A^ 
amounts to a real value for physical observables. It is 
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understood in view of the relation, 



det[ 7M (d M - ig(A fl ) c ) + j Af i q + m q ] 

= |det[7 M D A ' + 74Mg + ■ 



(4) 



We see clearly that the real (imaginary) part of the deter- 
minant is C-even (C-odd). For a C-even (C-odd) observ- 
able, thus, the imaginary (real) part of the determinant 
vanishes after integration over A^ 1 . Accordingly the gen- 
uine problem stems from that the real or imaginary part 
of the determinant may change its sign depending on the 
configuration A**. 



III. MODEL STUDY 

We will analyze a simple model to see the sign problem 
occurring in the mean-field level. As a practice to study 
the model, we will make the mean-field approximation 
in the SU(2) case for which we do not have to face the 
sign problem. Wc will then observe the sign problem in 
the SU(3) calculation and attempt the phase reweighting 
method to deal with the complex phase of the fermion 
determinant. 



A. Dense-Heavy Model 



the hopping parameter 1/ '(2m q a) and the gauge invariant 
chemical potential factor e fiqa [30(. After N T hops, a 
quark winds around the temporal circle. It results in 
the weight e for the quark excitation represented by the 
Polyakov loop L in the fundamental representation, 



= n u 4 (x,x 4 ) 

r l/T 

Pexp 



ig I dx 4 A^(x, Xi) 
o 



(7) 



The first line is the expression of the Polyakov loop in 
terms of link variables on the lattice and second is in 
the continuum. It does not matter whichever expression 
we use since the SU(JV C ) matrix L plays the role of the 
dynamical variable in our model and we will not return 
to its definition. 

The calculation of the determinant in color space is 
explicitly doable. After all we have the fermion action, 



-S f [L] _ 



■2e£] 



(8) 



in the color N c — 2 case (and Nf = 4 is implicit as we 
already noted) and 



Sf[L] _ 



= Y[[l + e 3 + 3e£ + 3e 2 t] 



(9) 



We will closely analyze a lattice model with dense 
heavy quarks [26|. Let us consider QCD in the limit 
of [i q — > oo so that we can drop antiquarks. We shall 
simultaneously take another limit of m q — > oo which ren- 
ders all quarks static. Under such limits we can evaluate 
the staggered fermion determinant exactly to reach the 
fermion action, 



-St[L] 



>[det(l + ei)] 



Nf/4 



(5) 



where Nf is the number of flavors. Not to go into subtlety 
of the flavor counting inherent to the staggered formal- 
ism, which is not of our interest, we set Nf = 4 through- 
out this paper. 

We take those two limits in a way characterized by the 
parameter e ranging from zero to infinity; 



\2m a a/ 



(6) 



We will often call this model parameter as the "density" 
parameter because e has strong correlation to the quark 
number density as seen in Fig. [1] for the SU(2) case and 
Fig. 0for the SU(3) case. Here N T is the number of the 
lattice sites in the temporal direction, that is, the inverse 
temperature 1/T is given by N T a with the lattice spacing 
a. In this model quarks are allowed to propagate only 
in the positive temporal direction. Each time a quark 
with mass m q travels by one temporal lattice, it picks up 



in the color N c = 3 case. In the above expressions we 
employed the traced Polyakov loop defined as 



1 



-trL, 



(10) 



where tr is taken in fundamental color space. We remark 
that the C-transformation changes I to £* and vice versa, 
where I is real for N c — 2 and generally complex for 
A c + 2. 

It is obvious that the A c = 2 determinant © is real 
and positive semidefinite for any e because 1 + e 2 + 2e£ > 
(1 — |e|) 2 > 0, while the A c = 3 expression © suffers 
the sign problem for nonzero e; the determinant can be 
complex except when either e = (zero density), e = 1 
(half- filling), or e — > oo (full-filling). We can intuitively 
understand why the SU(3) determinant becomes real for 
e = 1. One-quark excitation represented by £ and two- 
quark excitation that is equivalent to one-antiquark exci- 
tation represented by £* occur with the common weight 
e = e 2 = 1 because of half-filling (see n = 0.5 at £ = 1 
in Figs. [1] and [S]). Therefore, in effect, the system has 
equality in number of quarks and (effective) antiquarks 
just like in the zero density case. It is not an escape from 
the sign problem, however. In the A c = 3 case £ + £* can 
take a value ranging from —1 to 2, and thus the determi- 
nant at e = 1, i.e. det(l + L) oc 2 + 3(1 + £*) is real but 
can be negative for -1 < £ + £* < -2/3. 

We note here one more important feature of the model. 
The fermion determinant is invariant under the duality 
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transformation [2i 



1/e and £ 



(11) 



by which it is sufficient for us to investigate the model in 
the region e 6 [0, 1] and the outer region e <E (1, oo] can 
be deduced by means of the duality. 

Regarding the gluodynamics, we assume the nearest 
neighbor interaction between the Polyakov loops, 



s e {L} = ~NijJ2me*(y)- 



(12) 



This action looks pretty simple and still reproduces the 
fundamental nature of the phase transition in the pure 
gluonic sector; the action (fT2"l) leads to a second-order 
phase transition for the SU(2) case and a first-order phase 
transition for the SU(3) case [3l[, which is in agreement 
with the lattice QCD simulation and the theoretical ex- 
pectation from center symmetry 

One can interpret J as a model parameter specify- 
ing the "temperature" of the system. In the strong 
coupling expansion, in fact, J is related to T through 
J = exp[— aa/T] where a is the string tension. In this 
work we shall leave J as a model parameter as it is, for 
we do not want to introduce any further modeling into 
our analyses. In other words, our aim is to form a model 
not to imitate QCD itself but to mimic the QCD sign 
problem. 

The effective action that defines our model is eventu- 
ally given by 



S[L]=S g [L]+S { [L] 



(13) 



with the "density" parameter e contained in Sf[L] and 
the "temperature" parameter J in S g [L]. 



B. Mean Field Approximation 

In the finite-temperature field theory the free energy is 
evaluated in the functional integral form like the effective 
action as 



-f-V/T 



VLe 



-S[L] 



(14) 



We make use of the mean-field technique to approximate 
the free energy. Our ansatz for the mean-field action 

is EES 



S mf [L] = -^[^)+r(z)] 



(15) 



In case of pure gluonic theories, the mean-field x is sim- 
ply proportional to the Polyakov loop expectation value; 
x = 12N 2 J(£). In the presence of fermionic contribu- 
tions, however, there is no simple relation between them. 
Besides, (£) and (£*) in the SU(3) case have different 
dependence on [i q . One might think that there should 



be two independent mean-fields to deal with differing (£) 
and {£*). This idea has much to do with the sign problem 
actually, and we will come back to this point later. 
Then the mean-field free energy can be estimated as 



f mf (x)-V/T = (5[L]-5 mf [L]) mf -ln J VL, 



(16) 



where the average (• • • } m f is taken by the mean-field ac- 
tion S m {[L]. Roughly speaking, the first part corresponds 
to the internal energy and the logarithmic part is the en- 
tropy. We fix x so as to minimize / m f(x). Once S m i[L] is 
known with x determined, the expectation value of any 
observable 0[L] as a function of the Polyakov loop can 
be estimated by the group integration over L with the 
mean-field action S m i[L]; 

[ dLO{L}e- s ™ l[L] 
(0[L}) ~ (0[L]) mi ee J— . (17) 



For instance, the quark number density per color degrees 
of freedom is available by calculating 



_L df ~ £ / dSi 



N C V\ de 



(18) 



mf 



One can calculate the Polyakov loop susceptibility \ m 
the same way which reflects information of the deconfinc- 
ment phase transition. If one directly uses 0[L] — I 2 in 
the mean-field approximation, however, nothing becomes 
singular at the critical point unless the fluctuation of I 
is taken into account. Equivalently one can estimate the 
susceptibility from the inverse curvature of the effective 
potential because the one-loop fluctuation leads to the 
trace of propagator for (£ 2 ), that is, the inverse screening 
mass. Since (£) and (£*) are uniquely determined given x 
is fixed by the free energy (|16p . one can regard f m {(x) as 
a function of (£) or (£* } . The Polyakov loop susceptibility 
is then 



(d 2 f, 



mi 



dx 



d 2 f m 
dx 2 



(19) 



where we used df m f/dx = 0. Of course, the suscepti- 
bility defined in terms of £* is available with (d(£)/dx) 2 
replaced by (d{£*)/dx) 2 . The difference thus lies only in 
the nonsingular coefficient we are not interested in. For 
our purpose it is rather convenient to focus on the sin- 
gular part alone discarding the difference of (£) and (£*). 
Hence, we define the susceptibility as 



X 



d 2 x 



(20) 



for presenting our numerical results. 

Now that we finish explaining our approximations and 
computational procedures, let us step forward to the 
model analysis. 
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Density Parameter e Temperature Parameter / 



FIG. 1: Correlation between the density parameter e and the 
quark number density per lattice site (i.e. density in the unit 
of the lattice spacing) divided by N c = 2. It is obvious from 
the figure that e = is the zero-density (n = 0) state and 
e = 1 is the half-filling (n = 0.5) one. 



C. SU(2) Results 



FIG. 2: SU(2) (traced) Polyakov loop (I) as a function of the 
coupling J at various density parameters; e =0 (solid), 0.1 
(dotted), 0.2 (short-dashed), 0.5 (dashed), and 1.0 (dotted- 
dashed). The second-order phase transition at e = occurs 
at J c ~ 0.083. The transitional behavior is gradually smeared 
by the center symmetry breaking terms in the fermion deter- 
minant as e grows larger. 



We consider the model in the color SU(2) case first, as 
we mentioned, to see how the mean-field approximation 
works apart from the sign problem. 

Figure [1] shows the results for the quark number den- 
sity n as a function of the density parameter e using 
Eq. (fT8|) . For e > 1 the duality relation n(e) = 1 — n(l/e) 
enables us to deduce the number density. We can im- 
mediately confirm from this plot and the duality relation 
that the density parameter specifies the quark number 
density uniquely which monotonously approaches unity 
as e — > oo. It should be noted that the half-filling n = 0.5 
realizes at e = 1 and there is no J dependence at all then. 

Let us look at the phase transition seen in the Polyakov 
loop behavior with increasing J. The deconfinement 
phase transition in the SU(2) pure gluonic theory is 
known to be second-order bel onging to the same univer- 
sality class as the Ising model jl3l. [TS[32l |. In our model at 
e = we have a continuous transition at J = J c ~ 0.083 
as indicated by the solid curve in Fig. [5] The presence of 
dynamical quarks acts on the Polyakov loop variable as 
an external field breaking center symmetry. In fact, the 
results at nonzero e in Fig. [5] are not of transition but of 
crossover. 

We plot the density dependence of the Polyakov loop 
behavior in Fig. [3] The density effects generally tend 
to make the Polyakov loop larger, and eventually, the 
Polyakov loop becomes insensitive to the density in the 
large J (i.e. high temperature) region. It is because both 
the temperature and the density break center symme- 
try spontaneously and explicitly, respectively, having the 



Polyakov loop saturated. Therefore, the J dependence is 
less for larger e and the e dependence is less for larger J. 

Our mean-field outputs are to be compared with the 
lattice simulation in Ref. [26| ; our FigsJT] and [3] corre- 
spond to Figs. 1 and 2 presented in Ref. [26[ , respectively. 
We cannot expect a quantitative coincidence because our 
ansatz for the pure gluonic action S S [L] is only a crude 
approximation of QCD and besides we neglect the renor- 
malization of the Polyakov loop in the mean-field treat- 
ment. Nevertheless, the agreement turns out to be sur- 
prisingly good beyond our expectation if we treat the 
model parameter J as a fitting parameter incorporating 
the undetermined effect of the Polyakov loop renormal- 
ization, which is implied by the ansatz (|12p . In such a 
way, we can fix J = 0.0042 and J — 0.04 to reproduce 
the SU(2) Polyakov loop only at e = 1 for A/g 2 = 2.0 
and 4/g 2 = 1.5, respectively. We would emphasize that 
we did not use the data of the Polyakov loop at e ^ 1 
and not the data of the number density at all. Never- 
theless, as clearly seen from the comparisons in Figs. [4] 
and our numerical results fit all of the lattice data 
pretty well. We can conclude from this observation that 
the main QCD corrections to our ansatz (fT2]) could be 
absorbed into the renormalization of the coupling alone. 
We are now confident that the mean-field treatment is a 
fairly acceptable approximation for this type of problem. 

Finally let us check that the susceptibility ([20]) diverges 
at e = and J = J c . Figure [5] shows the susceptibility as 
a function of J at various e. In the plot x becomes greater 
with increasing J because we did not include d{t)/dx in 
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FIG. 3: SU(2) Polyakov loop {£) as a function of the density 
parameter e at various temperature parameters; J—0 (solid), 
0.05 (short-dashed), 0.1 (dashed), and 0.2 (dotted-dashed). 
The Polyakov loop becomes insensitive to e as J goes larger, 
which is consistent with Fig. [5] in which the results at various 
e converge at large J. 




Density Parameter e 

FIG. 4: Comparison to the number density measured on the 
lattice at 4/g 2 = 2. ( circle) and 4/ g 2 = 1.5 (square) taken 
from Fig. 1 in Ref. 0. 



our definition of \ in Eq. (I20[) . It is intriguing to remark 
that the e = 0.1 result is not really critical in view of Xi 
while the crossover at e = 0.1 in Fig. [5] looks rather close 
to a phase transition. Actually x is a more informative 
quantity to judge how critical the crossover is in fact. 

In summary of the SU(2) Polyakov loop dynamics at 
finite temperature and density, we shall depict a three- 
dimensional plot of (£) in Fig. as a function of the 



FIG. 5: Comparison to the SU(2) Polyakov loop measured 
on the lattice at 4/gr 2 = 2.0 (circle) and 4/<? 2 = 1.5 (square) 
taken from Fig. 2 in Ref. [2f|. 
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Temperature Parameter / 

FIG. 6: Susceptibility relevant to the SU(2) Polyakov loop I as 
a function of J at various density parameters; e =0.0 (solid), 
0.1 (short-dashed), 0.2 (dashed), and 1.0 (dotted-dashed). 
The e = result has a divergence at J = J c characteristic 
to the second-order phase transition. 



"temperature" J and the "density" e. We immediately 
see general tendency that (£) grows up with increasing J 
and e. 

Before closing this subsection we will comment a bit 
on the adjoint Polyakov loop whose definition is 



£ 



1 



adj 



-tri adj 



i 



(21) 



7V C 2 - 1 iV2 

The adjoint Polyakov loop is no longer an order param- 
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Temperature Parameter / 

FIG. 8: Dependence of the SU(2) adjoint Polyakov loop 
(^adj ) on the temperature parameter J at various density pa- 
rameters; e=0.0 (solid), 0.1 (dotted), 0.2 (short-dashed), 0.5 
(dashed), and 1.0 (dotted-dashed). 



eter because the adjoint representation is not faithful by 
the center of the gauge group. This quantity has been 
closely argued in Ref. [33] and there still remain sub- 
tleties. Not in the SU(2) but in the SU(3) case, it has 
been found that (£ a dj ) takes an infinitesimally small value 
in the confined phase. One possible explanation would be 
that the group integration over L is important in the cen- 
ter symmetric regime where nonperturbative phenomena 
like color confinement are relevant (see also Ref. 34]). 
As a matter of fact, the SU(-/V C ) group integration over 



£ a dj turns out to be zero, if we are allowed to disregard 
the dynamics. Even with the dynamics taken into ac- 
count within the mean-field approximation, as shown in 
Fig. [HJ the argument should hold so (^ a dj ) ~ in the low 
temperature side. However, the adjoint Polyakov loop 
behavior in our study should be understood only up to 
a qualitative level. It is pointed out in Ref. [33j that 
the renormalization for ^ a( jj is significant, which is not 
considered in our present treatment. 



D. SU(3) case 

In case of iV c = 3 we have to tackle the sign problem. 
We are not capable of solving QCD exactly like the lattice 
simulation, so one might think that within the framework 
of approximations one is allowed to impose the mean- field 
ansatz (|15p to get some results anyhow. The free energy 
after the integration over L could be real in terms of x. It 
seems to work at least as a rough estimate that is worth 
trying first. 

The serious flaw in such a simple strategy is that (£) = 
(£*) is inevitably concluded. It is, however, contradict to 
the lattice results and the model analyses [25[ where 
(I) ^ {£*) has been observed at finite density. If the 
mean-field ansatz (|15[) is extended to having two variables 
x and y in order to take account of the difference between 
(£) and (£*), the price to pay is that the mean- field free 
energy is not convex. We will revisit this issue latter. In 
any case, though the appearance might be unalike, the 
difficulty of the sign problem is conserved even in the 
mean- field approximation unless the difference (£) ^ (£*) 
is neglected. 

In our work we shall elucidate that the "phase 
reweighting method" is one way to resolve these diffi- 
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culties. We should, however, note that the reweight- 
ing method in the present context is one approxima- 
tion scheme unlike in the lattice simulation. That is, 
the reweighting method is expected to be precise if the 
number of configurations is infinitely large, and thus the 
lattice simulation with infinite number of configurations 
generated could provide us with the exact answer in prin- 
ciple, while the mean-field approximation picking up only 
the most dominant configuration cannot. 



1. Method 

The point of the method is that we decompose the 
fermion determinant into one part that gives the pos- 
itive semidefinite probability and the other part that is 
regarded as the observable whose average is taken by con- 
figurations. 

The complex fermion determinant consists of the C- 
even magnitude and the C-odd phase. Accordingly, the 
fermion action can be rewritten as 



S i [L) = Sr s [L]+i®[L] 

where 

S f mag [L] = - ^ ln| 1 + e 3 + 3e£ + ie 2 t \ 

X 

G[L] = -^2 ai 'g(! + fi3 + 3 <^ + 3e2 ^ 



(22) 

(23) 
(24) 



With these definitions we approximate the expectation 
value of 0[L] by the one obtained as follows; 



(0[L]) 



(0[L], 



-i&[L] 



mf 



(25) 



mf 



Here S m f[L] or x is fixed from the free energy with the 
action S e [L] + S™ as [L], so that x encompasses the infor- 
mation of S'™ as [L] implicitly. This scheme is the same 
as what has been adopted in the lattice simulation in 
Ref. J2|. 

Here we would draw attention to a related work; in 
Ref. [HI the correlation between lm£ and Q[L] was in- 
vestigated in the lattice simulation, which is apparent in 
our case from the expression (j2"4")l . 



2. Results 

Let us start with checking the monotonous correlation 
between e and n as in the SU(2) case. This allows us to 
regard increase (or decrease) in the density parameter e 
as increase (or decrease) in the quark number density n. 
The positive correlation is obvious from the results we 
show in Fig. [5] The "temperature" or J dependence is 
slightly greater than the SU(2) results. We can give a 
possible account for this as follows; in the confined phase 
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FIG. 9: Correlation between the density parameter e and the 
quark number density n per lattice site divided by iV c = 3. 
The gross feature is similar to Fig. [1] e = is the zero density 
state (n = 0) and e = 1 is the half-filling state (n = 0.5) in 
this SU(3) case as well as in the SU(2) case. 




Temperature Parameter / 

FIG. 10: SU(3) (traced) Polyakov loop (I) as a function of 
the temperature parameter J at various density parameters; 
e =0 (solid), 0.1 (dotted), 0.2 (short-dashed), 0.5 (dashed), 
and 1.0 (dotted-dashed). The first-order phase transition at 
e — and e = 0.1 occurs at J c == 0.132 and J c = 0.123, 
respectively. 



at small J, the quark number density is suppressed as 
compared with high J results. This suppression comes 
from the group integration that forces the thermally ex- 
cited particles to be not quarks but (nearly) color-singlet 
baryons consisting of -/V c quarks [Lsl . [3t| . In general larger 
N c leads to stronger suppression by heavier excitation 



0.5 1 

Density Parameter e 

FIG. 11: SU(3) Polyakov loop {£} as a function the density 
parameter e at various temperature parameters; J—0 (solid), 
0.05 (dotted), 0.10 (short-dashed), 0.13 (dashed), and 0.2 
(dotted-dashed). The gross feature is similar to the SU(2) 
case in Fig. [3] 



quanta. Therefore, the stronger J dependence presented 
in Fig. [^originates from the stronger suppression at small 
J. The suppression is physically interpreted as effective 
tendency toward confinement (2fj| . 

Figure [10] is the Polyakov loop as a function of J cor- 
responding to the SU(2) result in Fig. [5] and to be com- 
pare d q ualitatively with the lattice result of Fig. 7 in 
Ref. [26[ . We find a first-order phase transition for e = 
at J = J c = 0.132 and for e = 0.1 at J = J c = 0.123. 
The effect of nonzero e smears the transitional behav- 
ior and the phase transition eventually ceases to be of 
first-order at a certain e. The end-point of the first-order 
phase boundary is a second-order critical point called the 
critical end-point, which is of much interest in attempts 
to clarify the QCD phase diagram [37], HI]. We have 
crossover at larger e. The global picture is well consis- 
tent with what has been already clarified in the Potts 
system as a toy model of finite temperature and density 

QCD [nm. 

We plot the "density" dependence of the SU(3) 
Polyakov loop in Fig. [TT1 which is the SU(3) counterpart 
of Fig. [3] The SU(2) and SU(3) results are qualitatively 
similar except for that the Polyakov loop is suppressed 
at small J and e just as we found in the quark number 
density. It would be interesting if we could compare our 
results with the lattice data, but unfortunately, the SU(3) 
data as a function of e is not available from Ref. [26| . We 
cannot argue the J (or 6/g 2 ) dependence because it in- 
volves unknown renormalization effects. 

In the phase reweighting calculation we can see how 
(£) and (£*) become distinct at fx q =^ 0. The observable 
£ — £* is C-odd and so the imaginary part of the fermion 




Temperature Parameter J 

FIG. 12: Difference of the SU(3) Polyakov loops {£) and {£*) 
as a function of the temperature parameter J at various den- 
sity parameters; e =0 (solid) which is zero entirely, 0.1 (dot- 
ted), 0.2 (short-dashed), and 0.5 (dashed). 
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FIG. 13: Difference of the SU(3) Polyakov loops (£) and 
{£*) as a function of the density parameter e at various tem- 
perature parameters; J =0 (solid), 0.05 (short-dashed), 0.1 
(dashed), and 0.2 (dotted-dashed). The thin curves represent 
the results from the Taylor expansion method in the confine 
phase at J < J c and in the deconfined phase at J = 0.20 > J c 
which is almost overlaid on the result from the phase reweight- 
ing. 



determinant is responsible for a nonvanishing difference. 
When e is small in the fermionic determinant © the 
imaginary part comes from lme£ oc e(£ — £*). Conse- 
quently the expectation value of the difference is propor- 
tional to e((Im£) 2 )o where (• • • )o is taken at zero den- 
sity [25[ | . In Fig. [T^] we present our numerical results for 
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the difference {£) — {£*) as a function of J. The difference 
is trivially zero at e = and e — 1 where the fermion de- 
terminant is real. As long as the density parameter stays 
smaller than e ~ 0.5, a larger density parameter e leads 
to a bigger difference. For example, we find the differ- 
ence at e = 0.5 as large as (£) — (£*) = —0.076 which is 
comparable to (£) = 0.073. 

One can intuitively understand why (£*) is greater than 
(£) at nonzero fi q , which agrees with what has been ob- 
served in the lattice simulation [n}. It is because, as 
discussed also in Ref. [lljjj the presence of quarks at fi- 
nite density enhances the screening effect for antiquaries 
so that the antiquark excitation costs less energy. 

The dense-heavy model, or e~ St , originally takes a 
form of power series in e as seen in the expression (|9l). 
We have thus performed the Taylor expansion also to 
estimate the Polyakov loop difference. The expectation 
value of the coefficient of the e series is calculated at zero 
density e = 0. In this static model the fermionic contri- 
bution is just vanishing at e = 0, and so the Taylor ex- 
pansion method gives rise to the identical output for any 
J < J c in the confined phase where x remains zero. In 
the deconfined phase J > J c the J dependence is brought 
in by nonzero x. Figure [T3l shows the difference (£) — {£*) 
as a function of e at various J with the Taylor expansion 
results in the confined and deconfined phase. We can see 
excellent agreement between two methods from the com- 
parison at J = 0.00 and J = 0.20. However, the Taylor 
expansion cannot reproduce the results at < J < J c at 
all. The lesson from our model study is that the Taylor 
expansion in terms of density breaks down in the pres- 
ence of the first-order phase transition with respect to 
temperature. Still, in the deconfined phase at high tem- 
perature, the expansion is validated. We do not think 
that this finding is trivial in our model. In the realistic 




£=1.0 



Q| , . , , , , , , , 1 

0.1 0.2 

Temperature Parameter / 

FIG. 15: Susceptibility relevant to the SU(3) Polyakov loop 
t as a function of J at various density parameters; e = 0.0 
(solid), 0.1 (short-dashed), 0.2 (dashed), and 1.0 (dotted- 
dashed) . 



QCD lattice simulation with 2+1 flavors, the zero den- 
sity result is most likely crossover, and thus the Taylor 
expansion method should be reliable at all temperatures. 

As we showed in the previous subsection, we shall 
present the three-dimensional plot of the SU(3) Polyakov 
loop as a function of the temperature parameter J and 
the density parameter e in Fig. [Mj The figure is quali- 
tatively consistent with the lattice result shown in Fig. 5 
in Ref. [||. 

We shall turn to the susceptibility \ to examine the 
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Temperature Parameter / Density Parameter e 



FIG. 16: SU(3) adjoint Polyakov loop (4idj) as a function of 
the temperature parameter J at various density parameters; 
e=0.0 (solid), 0.1 (dotted), 0.2 (short-dashed), 0.5 (dashed), 
and 1.0 (dotted-dashed). 



FIG. 17: Phase of the fermion determinant Ke(e~'°) per lat- 
tice site as a function of the density parameter e at various 
temperature parameters; J =0 (solid), 0.05 (short-dashed), 
0.10 (dashed), and 0.20 (dotted-dashed). 



phase transition more closely. Because the phase transi- 
tion at small e is of first-order, \ jumps discontinuously 
at J — — J, c as seen in Fig. [J3J The susceptibility \ grows 
as the model parameters approach the critical end-point 
where the second-order phase transition takes place. Fig- 
ure [m implies that the zero density (e = 0) result is sig- 
nificantly affected by the critical end-point which should 
be located nearby at small e. In view of Figs. IBI and [Tol 
the critical region in the SU(3) case is wider than the 
SU(2) case. 

Let us comment on the adjoint Polyakov loop here 
again. As we mentioned in discussions on the SU(2) re- 
sults, the adjoint Polyakov loop should require proper 
renormalization beyond the mean- field treatment. We 
show the adjoint Polyakov loop result in Fig. [TB] just be- 
cause we can do that. The first-order phase transition 
is located at the same point; J = J c ~ 0.132, of course. 
The group integration over L suppresses at low J 

even for large e, which makes the crossover in Fig.lTfillook 
shaper than in case of the fundamental Polyakov loop in 
Fig. Ql 

Finally we present the results for the expectation value 
of the phase factor of the fermionic determinant, e~ lB . 
We plot Re(e~* e ) as a function of e in Fig. [T71 where 
6 is the phase at each lattice site; 9 = — arg(l + e 3 + 
3e£ + 3e 2 l*) (i.e. 6 = Ex 6 *-) Comparing it with Fig. 9 
in Ref . [26[ , we can surely check that our results qualita- 
tively reproduce the lattice data. For more quantitative 
arguments, let us put the volume 6 3 = 216 of the lattice 
simulation in Ref. [26[ into our results. The expected 
phase factor (e~ tB ) can be, as a rough estimate, approx- 
imated as ((e~ l9 )) 216 . For instance, our J = result has 
a minimum at e = 0.61 where (e.~ l9 ) — 0.977, and we 



get 0.977 216 = 0.0066. The minimum value in Fig. 9 in 
Ref. is read as of order 0.01, which is not quite far 
from our estimate viewed on the logarithmic plot. 



E. Mean-Field Free Energy 

Here we will pursue another strategy to deal with the 
difference between (£) and (£*) as a double-check. It is 
possible to extend the mean-field ansatz as 



Smf[i] = -f I)[W(2)]-5 ^2[l(S)-t(x)] , (26) 



and then we can compute the mean-field free energy from 
Eq. (|16p as a function of x and y. After the integration 
over L, the free energy f m f(x, y) is a real function of real 
variables x and y. We can thus fix the mean-fields by 



Of, 



mf 



dx 



Of 



mf 



(x,y) = (x ,yo) 



dy 



0. 



(27) 



(x,y)=(x ,yo) 



Nonzero yo appears in the presence of nonzero fi q . It 
turns out that around (xo,j/o) the free energy f m f{x,y) 
has a minimum in the x direction, while it has a maxi- 
mum in the y direction. That is, the solution to Eq. (|27|) 
is a saddle-point of f m f(x,y), which is consistent with 
Ref. . The instability with resp ect to y should be 
a remnant of the sign problem [4l| We here point out 
that this instability also exists in the zero density limit; 
fmi{x, y) has a saddle-point at y = even at zero den- 
sity. It should be instructive to take a closer look at how 
the saddle-point arises even at zero density. If we expand 
the free energy in terms of (£) and (£*), the leading term 



12 



o 



-0.05 



£■=0.1 
£•=0.5 
£•=0.1 (MF) 
£•=0.5 (MF) 




0.1 0.2 

Temperature Parameter J 

FIG. 18: Comparison to the difference {£) — (I*) estimated by 
the saddle-point of the mean-field free energy at various den- 
sities; e =0.1 (solid and dashed) and 0.5 (dotted and dotted- 
dashed) . 



dependent on {£) and {£*) is quadratic ~ (£)(£*) as explic- 
itly seen in a simple estimate in Ref. • This form of the 
free energy implies an instability inducing {£} ^ (£*) be- 
cause it can written as (£)(£*) oc {{£)+{£* )) 2 -{{£)-{£* )) 2 . 

The situation is somewhat analogous to thermodynam- 
ics in the finite density NJL-model calculation. If the free 
energy / is calculated as a function of the renormalizcd 
chemical potential \i q [38| . then the value of \i q is fixed 
by the condition df/d/j, q = which corresponds to not 
a minimum but a maximum of / as a function of \i q . 
This is not problematic, however, because the condition 
df /dj-iq — means a constraint equation of number den- 
sity. Likewise, we might as well think that the determina- 
tion of y in the Polyakov loop dynamics is not energetic 
but the second equation of ([27]) has something to do with 
a constraint equation of number density of gauge charge, 
namely, the Gauss law constraint. If this conjecture is 
the case, though the rigorous proof is beyond our current 
scope, the saddle-point nature of the free energy is no 
longer an obstacle to fix x and y. 

Let us see what comes out if we calculate the Polyakov 
loop using the ansatz (|2l)|) with the solution to Eq. (|2"T]) 
assuming that the second equation of (f27|) stems from a 
constraint associated with gauge dynamics. In Fig.fTSlwe 
show the difference (£) — (£*) as a function of J. The solid 
and dotted curves at e = 0.1 and e = 0.5 respectively are 
just the same as already presented in Fig.[l2j The dashed 
and dotted-dashed curves are the counterparts derived 
from the mean-field energy with x and y. 

In view of two e = 0.5 curves, on the one hand, 
the mean-field result is entirely consistent with what we 
found in the phase reweighting method. On the other 
hand, the coincidence is not very good for the e = 0.1 
results except for J ~ and J > J c . This discrep- 



ancy comes from a singularity of the dense-heavy model 
located at e = 0; when e is small, the fcrmion action 
in the dense- heavy model is approximated as Sf[L] ~ 
— 3e £{x). Hence, nonzero e tends to align the vacuum 
into the direction of x — y (see Eq. (|26)) ), and the model 
does not reduce to a pure gluonic theory with y = 
smoothly in the limit of e — > + . This situation makes a 
sharp contrast to finite density QCD. In the strong cou- 
pling expansion, for a simple example, the fermion action 
is Sf[L] ~ H Ea?[A£) e/ia/T + Z*{x)e-^/ T ] . In the limit 
of fi q — > 0, therefore, the fermionic action acts as an ex- 
ternal field toward x ^ and y — 0, so that the vacuum 
alignment of a pure gluonic theory is smoothly retrieved 
in the /j, q — > and H — > + limit. 

Apart from the discrepancy inherent to the singular 
behavior of the dense-heavy model near e = 0, the mean- 
field free energy leads to the results in accord with the 
phase reweighting method. This is an indirect evidence 
for that the saddle-point is not harmful actually, as we 
conjectured. 

We shall comment on a possible clue to resolve the 
sign problem based on what we have seen here. We 
changed the sign problem into a form of the saddle-point 
of the mean-field free energy. The saddle-point appears 
harmless from our analyses, presumably stabilized by the 
Gauss law, and then the sign problem is resolved. Of 
course, it is highly nontrivial how to map this analyti- 
cal procedure to the lattice QCD simulation. Still, we 
would point out that the clustering method developed 
in Ref. [39| is a philosophically similar idea along this 
line; partial integration over the cluster domains wipes 
away the sign problem, just like seen in our mean-field 
free energy free from the sign problem after the group 
integration. 

The field-theoretical approach to reveal the relation 
between the saddle-point of the free energy and the Gauss 
law is beyond our current scope, but it definitely deserves 
further investigation. 



IV. SUMMARY 

We investigated the dense-heavy model and observed 
the sign problem at finite density within the framework of 
the mean-field approximation. We calculated the quark 
number density, the fundamental and adjoint Polyakov 
loop, the susceptibility, and the phase of the fermion de- 
terminant as a function of the model parameters specify- 
ing the temperature and the density. All the mean-field 
results are reasonable and even in quantitative agreement 
with the lattice data. 

In the environment free from the sign problem in the 
SU(2) case we found that the mean-field approximation 
goes better than expected. Our mean-field results nicely 
reproduced the quark number density and the Polyakov 
loop measured on the lattice once the unknown Polyakov 
loop renormalization is fixed by fitting. Then, we pro- 
ceeded into the SU(3) case where the sign problem is 
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relevant. 

We saw that the approximation scheme which is capa- 
ble of describing the different fx q dependence of (£) and 
(£*} cannot avoid the sign problem even at the mean- 
field level. We applied the phase reweighting method 
as a practical resolution in order to handle the complex 
fermion determinant. We acquired (£) — {£*} as a function 
of the density parameter or the temperature parameter. 
So far, there are not many lattice data available for this 
quantity, but our results (£*) > {£) are consistent with 
other model studies as well as what has been obtained 
from the Taylor expansion method on the lattice. 

The important message from our work is the following. 
The sign problem may be a serious obstacle even in the 
mean-field model studies at finite temperature and den- 
sity. We would say, at least, that one should be careful 
enough when one deals with the Polyakov loop behav- 
ior at finite density. The chiral effective models with the 
Polyakov loop coupled are examples that definitely need 
more or less caution for application to the finite density 
problem. The phase reweighting method is one prescrip- 
tion in order to approximate the expectation value in 
a manageable way. This prescription is, however, only 



practical and not any solution to the sign problem. We 
would insist that one cannot solve the QCD sign problem 
until one can at least find a solution to the sign problem 
appearing even in that simple model setting. The con- 
verse is not necessarily true, though. 

In the future, we would like to apply other ideas than 
the phase reweighting method for the sign problem at 
the mean-field level. The bottom line here is, thus, that 
we have formulated an analytical testing ground to think 
about the sign problem, and we believe that our clarifi- 
cation would be useful for further developments. 
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